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Geomètric phases have stimulated researchers for its potential applications in many areas of 
science. One of them is fault-tolerant quantum computation. A preliminary requisite of quantum 
computation is the implementation of controlled lògic gates by controlled dynamics of qubits. In 
controlled dynamics, one qubit undergoes coherent evolution and acquires appropriate phase, de- 
pending on the state of other qubits. If the evolution is geomètric, then the phase acquired depend 
only on the geometry of the path executed, and is robust against certain types of errors. This phe- 
nomenon leads to an inherently fault-tolerant quantum computation. Here we suggest a technique 
of using non-adiabatic geomètric phase for quantum computation, using selective excitation. In a 
two-qubit system, we selectively evolve a suitable subsystem where the control qubit is in state 
|1), through a closed circuit. By this evolution, the target qubit gains a phase controlled by the 
state of the control qubit. Using these geomètric phase gates we demonstrate implementation of 
Deutsch-Jozsa algorithm and Grover's search algorithm in a two-qubit system. 

I. INTRODUCTION 



^ ■ Quantum states which differ only by a overall phase cannot be distinguished by measurements in 

, H : 

l 5t .quantum mechanics. Hence phases were thought to be unimportant until Berry made an important 
and interesting observation regarding the behavior of pure quantum systems in a slowly changing 
environment [1]. The adiabatic theorem makes sure that, if a system is initially in an eigenstate of the 
instantaneous Hamiltonian, it remains so. When the environment (more precisely, the Hamiltonian) 
returns to it's initial state after undergoing slow changes, the system acquires a measurable phase, 
apart from the well known dynamical phase, which is purely of geomètric origin [1] . Simón [2] showed 
this to be a consequence of parallel transport in a curved space appropriate to the quantum system. 
Berry's phase was reconsidered by Aharanov and Anandan, who shifted the emphasis from changes 
in the environment, to the motion of the pure quantum system itself and found that the for all 
the changes in the environment, the same geomètric phase is obtained which is uniquely associated 
with the motion of the pure quantum system and hence enabled them to generalize Berry's phase to 
non-adiabatic motions [3]. For a spin half particle subjected to a magnètic field B, the non adiabatic 
cyclic Aharanov-Anandan phase is just the sòlid angle determined by the path in the projective 



Hilbert space [3]. 

Yet another interesting discovery in the fundamentals of quantum physics was the observation that 
by accessing a large Hilbert space spanned by the linear combination of quantum states and by intel- 
ligently manipulating them, some of the problems intractable for classical computers can be solved 
efficiently [4,5]. This idea of quantum computation using coherent quantum mechanical systems has 
excited a number of research groups [6,7]. Various physical systems including nuclear magnètic reso- 
nance (NMR) are being examined to built a suitable physical device which would perform quantum 
information processing and quantum simulations [8-12]. Also, the quantum correlation inherently 
present in the entangled quantum states was found to be useful for quantum computation, commu- 
nication and cryptography [5]. Geomètric quantum computing is a way of manipulating quantum 
states using quantum gates based on geomètric phase shifts [13,14]. This approach is particularly 
useful because of the built-in fault tolerance, which arises due to the fact that geomètric phases de- 
pend only upon some global geomètric features and it is robust against certain errors and dephasing 
[13-16]. 

In nuclear magnètic resonance (NMR), the acquisition of geomètric phase by a spin was first 
verified by Pines et. al. [17] in adiabatic regime by subjecting a nuclear spin to an effective magnètic 
field which slowly sweeps a cone. A similar approach was adopted by Jones et. al. to demonstrate the 
construction of controlled phase shift gates in a two-qubit system using adiabatic geomètric phase 
[13,14]. Pines et. al. also studied the geomètric phase in non-adiabatic regime, namely the Aharanov- 
Anandan phase by NMR [18]. They used a system of two dipolar coupled identical proton spins 
which form a three level system. A two level subsystem was made to undergo a cyclic evolution in the 
Hilbert space by applying a time- dependent magnètic field, while geomètric phase was observed in the 
modulation of the coherence of the other two-level subsystem [18]. Recently, non-adiabatic geomètric 
phase has also been observed for mixed states by NMR using evolution in tilted Hamiltonian frame 
[19]. In the work reported here, we have adopted a scheme similar to that of Pines et. al. [18] to 
demonstrate construction of controlled phase shift gates in a two-qubit system using non-adiabatic 
geomètric phase by NMR. The scheme is easily scalable to higher qubit systems. The geomètric 
controlled phase gates were used to implement Deutsch-Jozsa (DJ) algorithm [20] and Grover's 
search algorithm [21] in the two-qubit system. To the best of our knowledge, this is the first 
implementation of quantum algorithms using geomètric phase. 



II. NON-ADIABATIC GEOMÈTRIC PHASE GATE 



Consider a two qubit system, which has four eigenstates 1 00) , 1 01) , 1 10) and |11). The two-state 
subsystem of 1 10) and 1 1 1) can be taken through a circuit enclosing a sòlid angle Q [18]. If the 
other dynamical phases are canceled during the process, these two states gain a non-adiabatic phase 
purely due to geomètric topology. Since the operation is done selectively with the states where first 
qubit is in state |1), this acts as a controlled phase gate where the second qubit gains a phase only 
when the first qubit is |1) [6,7,13] . 

The transport of the selected states through a closed circuit can be accomplished by selective 
excitation. Such selective excitation can be performed by pulses having a small bandwidth which 
excite a selected transition in the spectrum and leave the others unaffected [9,23-27]. In the following 
we consider geomètric phase acquired by two different paths in a Bloch sphere, respectively known 
as slice circuit and triangular circuit [18]. 

A. Geomètric phase acquired by a slice circuit 

In a slice circuit, the state vector cuts a slice out of the Bloch sphere, Figure l(a). The slice circuit 
can be achieved by two pulses A.B=(7r)[ ) 10 ^' 11 ^ (71)^^}^, where the pulses are applied from left 
to right. Here (tt)^ 10 ^' 11 ^ denotes a selective 7r-pulse on 1 10) <-> 1 1 1 ) transition with phase 6. The 
first (tt)^ 10 ^' 11 ^ pulse rotates the polarization vector of the subsystem through ir about the an axis 
with azimuthal angle 9 in the x-y plane (Fig. 1). The vector is brought back to its original position 
completing a closed circuit by the second 7r-pulse about the axis in the x-y plane with azimuthal 
angle (9 + rr + 0). The resulting path encloses a sòlid angle of 20, The operator of the two pulses 
can be calculated as, 

= exp{-i{I^\ ll) cos{e) + /l 10) ^ |11) sm(^))7r] 

exp[-i{I^ lV >cos{9 + 7T + 0) + I^"\ lV >sin{9 + tt + 0))vr] 
(l \ 
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where i"]. 10 )"! 11 ) and /I^H 11 } are the fictitious spin-1/2 operators [28] for the two-state subsystem of 
1 10) and |11), given by; 
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Note that the combined operator of the two pulses in Eq. [1] attributes a non-adiabatic geomètric 
phase proportional to the sòlid angle of the circuit traversed. However, the phase it attributed only 
to the two states where first qubit is in state |1). Since the selective excitation does not perturb 
the other transitions, the subsystem 1 00) <-> 1 01 ) do not gain any phase, This is analogous to the 
controlled phase gate where the second qubit acquires a phase controlled by the state of first qubit 
[7,13]. 

To demonstrate the operation of a controlled geomètric phase gate, we have taken the two qubit 
system of carbon-13 labeled chloroform ( 13 CHCl3), where the two nuclear spins 13 C and 1 H forms the 
two-qubit system. The sample of 13 CHC1 3 was dissolved in the solvent of CDC1 3 , and experiments 
were performed at room temperature at a magnètic field of B =11.2 Tesla. At this high-field the 
resonance frequency of proton is 500 MHz and that of carbon is 125 MHz. The indirect spin-spin 
coupling (the J-coupling) between the two qubits is 210 Hz. Starting from equilibrium, the 1 00) 
pseudopure state was prepared by spatial averaging method using the pulse sequence [19], 

1 



(n/3y x — G z - (tt/4); 



2J (*/4)-y-G z , 



(3) 



where the pulses were applied on the second qubit, denoted by 2 in superscript, which in our case 
is the proton spin. After creation of pps, a pseudo-Hadamard gate [29,30] was applied on the first 
qubit, which in our case was 13 C. The pseudo-Hadamard gate was implemented by a (7r/2)* where 
the superscript denotes the qubit and the subscript denotes the phase of the pulse [29,30]. This gate 



creates a an uniform superposition of the first qubit 1 00) + |10). The operation of the controlled 
phase gate would now transform the state into 1 00) + e^|10). For the slice circuit, the proton 
dynamical phase would vanish since the applied field is always orthogonal to the polarization vector, 
generating paral·lel transport [18]. However, the carbon coherence would undergo evolution due to 
the internal Hamiltonian during the pulses. Hence the pulse sequence of the gate was incorporated 
into a Hahn-echo [31,32] sequence of the form r — (Tr)^ — r, where the pulse sequence of the gate were 
applied during the second r period, as given in figure l(b). The intermediate (7r)-pulse refocuses 
inhomogeneity of the B field, the chemical shift of carbon and its J-coupling to the proton. However, 
to restore the state of the first qubit altered by the (7r)-pulse, the pulse sequence of Eq. [1] has to 
be supplemented by adding a {j)\ x pulse (figure l(b)), yielding the sequence: 
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(4) 



where the selective pulses were applied on the 1 00) <-> 1 1 ) transition to achieve the exact form of 
controlled phase gate. The selective excitation was obtained with Gaussian shaped pulses of 13.2 ms 
duration. The non-adiabatic geomètric phase was observed in the phase of 1 00) <-> 1 10) coherence. 
We have observed the geomètric phase for the slice circuit with various sòlid angles (20), each time 
varying the phase of the second selective (7r)-pulse. The corresponding spectra are given in figure 
2(c), where the 1 00) <-> |10) shows a phase change of é 1 ^. For = 0, there is no phase change and 
the peak is absorptive. With increase of 0, the phase of the peak changes and it becomes dispersive 
for = 7r/2, and subsequently, a negative absorptive for = ir. 

The three small lines in the spectra comes from the naturally abundant 13 C signal of CDCI3, which 
provide a reference. Since all dynamical phases due to evolution under chemical shift and J-couplings 
were refocused, the solvent 13 C signal is absorptive in all the spectra. However, solute 13 C signal 
gains phase because it is coupled to the protons, one of whose transition is taken through a closed 
circuit. This result thus provides a graphic display of geomètric phase by non-adiabatic evolution. 



To accurately read the phase angle of each spectrum in Fig. 2(c), a zero-order phase correction was 
applied to the spectra in Fig. 2(c), till the observed peak became absorptive. The change of phase 
of 1 00) <-> 1 10) coherence due to geomètric phase is plotted against the solid-angle (20), in figure 3. 
The graph in figure 3 shows the high fidelity of the experimental implementation of the slice circuit 
in this case. 



B. Geomètric phase acquired by a triangular circuit 



In the triangular circuit, the state vector traverses a triangular path on the Bloch sphere figure 
l(c) [18]. The sòlid angle enclosed by the triangular circuit of figure l(c) is <fi. The controlled phase 
shift gate can be implemented by the non-adiabatic phase acquired when the appropriate sub-system 
goes through this circuit. The pulse sequence for the circuit and the corresponding operator can be 
calculated, similar to that of the sliced circuit, as 
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The intermediate (0)l lo) ^ |n) pulse can be applied by the composite z-pulse sequence 



|io>Hii)(0)|io>Hii> 
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[33,34]. 



In the experiments, we have chosen 9 = 37r/2. The state of 1 00) + 1 10) was prepared and then the 
pulse sequence of figure l(d) was applied. Similar to the slice circuit, the sequence was incorporated 
in a Hahn-echo and the pulses were applied on the 1 00) <-> 1 1 ) transition. The operator of Eq.[5] 
transforms |00) + 1 10) to |00) + e"^/ 2 |10). The phase of the |00) <-> 1 10) was observed for various 
0, by changing the angle of the z-pulse and the phase of the last pulse in Eq.[5]. The spectra are 



given in figure 4. Once again, the peak changes from absorptive to dispersive and then to a negative 
absorptive in correspondence with the change of 0. 

However, there are two major differences between the spectra of figure 2(c) and 4(c). Note that 
after the phase gate, the state of the system is 1 00) + e t( ^\10) for slice circuit and 1 00) + e — *^/ 2 | 10) 
for triangular circuit. This is because the sòlid angle of the slice circuit if 20, whereas that of the 
triangular circuit is 0. Hence, in the slice circuit the coherences become a negative absorptive for 
= 7r, whereas in the triangle circuit the same observation is obtained for = 2n. Moreover, the 
phase of the pulses corresponding to the triangle circuit is chosen such that the sign of phase is 
opposite to that of the slice circuit. This difference is clearly reflected in the sign of the coherences 
between figure 2(c) and 4(c). A plot of the absolute value of observed phase change against sòlid 
angle is given in figure 5, whose high fidelity validate the use of such gates for quantum computing. 

III. DEUTSCH-JOZSA ALGORITHM 

Deutsch-Jozsa (DJ) algorithm provides a demonstration of the advantage of quantum superposi- 
tions over classical computing [20]. The DJ algorithm determines the type of an unknown function 
when it is either constant or balanced. In the simplest case, f(x) maps a single bit to a single bit. 
The function is called constant if f(x) is independent of x and it is balanced if f(x) is zero for one 
value of x and unity for the other value. For N qubit system, f(x±, X2, ...%) is constant if it is 
independent of Xi and balanced if it is zero for half the vàlues of Xi and unity for the other half. 
Classically it requires {2 N ~ l + 1) function calls to check if f{xi,X2,---x^) is constant or balanced. 
However the DJ algorithm would require only a single function call [20]. The Cleve version of DJ 
algorithm implemented by using a unitary transformation by the propagator Uf while adding an 
extra qubit, is given by [35], 

Uf 

\x 1 ,x 2 , ...x N )\x N+1 ) — > \x 1 ,x 2 , ...x N )\x N+1 © f(x 1 ,x 2 , ...x N )) (6) 

The four possible functions for the single-bit DJ algorithm are f 00 , fu, fio an d f m . foo(x) = for 
x =0 or 1, fn(x) = 1 for x =0 or 1, fio{x) =1 or corresponding to x =0 or 1, while foi(x) =0 
or 1 corresponding to x —0 or 1. The unitary transformations corresponding to the four possible 
propagators Uf are 
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(7) 



For higher qubits the functions are easy to evaluate using Eq. [6]. DJ-algorithm has been demon- 
strated using dynamic phase by several research groups [29,36-39]. 

The quantum circuit for single-bit Cleve version of DJ algorithm is given in figure 6(a) [36]. The 
algorithm starts with 1 00) pseudopure state. The pair of pseudo-Hadamard gates (ir / 2) y (ir / 2) 2 _ y 
create superposition of the form [(|0) + |l))/v / 2][(|0) — |l))/\/2]. Then the operator Uf is applied. 
When the function is constant, i.e. /(0) = /(l), the input qubit is in the state (|0) + \l))/y/2, else 
the function is balanced in which case it is in the state (|0) - \l))/y/2. Thus, the answer is stored in 
the relative phase between the two states of the input qubit. A final pair of pseudo-Hadamard gates 
(n / 2)\ y (ir / 2)y converts the superposition back into the eigenstates. The work qubit comes back to 
state |0), where as the input qubit becomes |0) or |1) corresponding to the function being constant 
or balanced. 

The operator of Uf 00 is identity matrix and corresponds to no operation. The operator of C// n 
can be achieved by a {ji) x pulse on the second qubit. In this experiment, unlike section II, we label 
proton as the first qubit and carbon as the second qubit, and consequently the (tt) x pulse was applied 
on the carbon. The U f 10 operator is a controlled-NOT gate which flips the second qubit when the 
first qubit is |1). This gate can be achieved by a controlled phase gate sandwiched between two 
pseudo-Hadamard gates on the second qubit [], Uf 10 — h — Cn(7r) — ftr 1 , where the controlled phase 
gate is of the form, 
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This precise form of controlled phase gate can be achieved by a recursive use of the phase gates 
demonstrated in section II. Since the gate A.B given in Eq.[l] attributes a phase e l<t> to the state 1 10) 
and to the state |11), we denote this gate as C w ((f)).Cn(— 0), where 



A.B = [C lo (0).C n ( 



/l \ 
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The phase gate Cn(0) can be constructed by a suitable combination of these gates, 
[C O o(-0/4).C , lo (0/4)] x [Coi(-0/4).C n (0/4)] x [C w (-<j>/2).C n (<f>/2)] 



(9) 
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Note that if performed in fault-tolerant manner by using non-adiabatic geomètric phase, the first 
gate requires a rotation of the transition 1 00) <-> 1 10} through a closed circuit. We have used the slice 
circuit, where it requires a sequence of two 7r-pulses, (^)é OO ^' 1O ^(^)é°+t^0 1 /4· Similarly, the second 
phase gate of Eq.[7] can be achieved by the pulse sequence (^e^^^i^g^^l^jl- Note that these 
two sequence is require pulsing of both the transitions of first qubit, 1 00) <-> 1 10} for the first gate and 
1 01} <-> 1 1 1 } for the second. Hence, they can be performed simultaneously by a couple spin-selective 
pulses (7r)e(7r)e+7r-0/45 where the pulses are applied on the first qubit (denoted by superscript). Thus, 



Cn(<l>) 



(vr)e.(7r)^_ 0/4 .(7r 



,110)^111) / \|10>«|11> 



(11) 



In this case = vr, and we have chosen 9 = 3n/2. The last two pulses are however transition selective 
pulses, which were incorporated into a refocusing sequence, r — (vr/2) 4 — r — (n/2) 2 — r — (7r/2) 4 — 
r — (tt/2)1, where the selective pulses were applied in the last r period, and the pulses were applied 
on the 1 00) <-> 1 1 } transition. It may be noted that the triangular circuit could have also used for 
the same purpose. The pseudo-Hadamard pulses on second qubit were achieved by h = (jt /2) 2 and 
h~ l = (n/2) 2 _ y pulses. 

The operator of U f 01 can be implemented in the similar manner by h — CqqÍjï) — h^ 1 , where C O o(0) 
can be implemented by 



c M ) = M^)í + , +4/4 .(4 10)HU U4Tt;£ 



(12) 



The equilibrium spectrum of the two qubits are given in figure 6(b). After creating the superposition 
from pps, applying the various Uf, and applying the last set of (ir/2) pulses, the spectra of proton 
and carbon were recorded in two different experiments by selective (vr/2) pulses after a gradient. The 
spectra corresponding to various functions are given in figure 6(c), (e), (g) and (i). The intensities 
of the peaks in the spectra provide a measure of the diagonal elements of the density matrix. The 
complete tomographed [40,41] density matrices in each case is given in figures 6(d), (f), (h) and (j). 
When Uf 00 and Uf n are implemented, the final state is 1 00) , and since the state of input qubit is 
|0), the corresponding functions foo and fu are inferred to be constant. Whereas in the case of Uf 01 
and Uf 10 , the final state of the system in |10). The state of input qubit being |1), the corresponding 
functions / i an d fio are balanced. Theoretically, it is expected that the density matrices will have 
only the populations corresponding to the final pure states. There were however errors due to r.f. 
inhomogeneity and relaxation. The deviation from the expected results are within 13%. 



IV. GROVER'S SEARCH ALGORITHM 

Grover's search algorithm can search an unsorted database of size N in 0{\fN) steps while a 
classical search would require O(N) steps [21]. Grover's search algorithm has been earlier demon- 
strated by several workers by NMR, all using dynamic phase [30,42-44,39]. The quantum circuit for 
implementing Grover's search algorithm on two qubit system is given in figure 7(a). The algorithm 
starts from a 1 00) pseudopure state. A uniform superposition of all states are created by the initial 
Hadamard gates (H). Then the sign of the searched state "x" is inverted by the oracle through the 
operator 

U x = I-2\x)(x\, (13) 

where U x is a controlled phase shift gate C x {ji). Cn(7r) and C 00 (n) gates were implemented by the 
pulse sequences given in Eq.[ll] and [12] respectively. The oracle for the other two states |01) and 
1 10) were implemented by the sequences, 

CM) = (14) 

where = ir, as required in our case. 
An inversion about mean is performed on all the states by a diffusion operator HU 00 H [21], where 

E/oo = J-2|00)<00|, (15) 



where U 00 is nothing but C 00 (7i), and was implemented by the pulse sequence of Eq.[12]. For an 
N-sized database the algorithm requires 0(y/~N) iterations of U X HU 00 H [21]. For a 2-qubit system 
with four states, only one iteration is required [30,42]. We have created a 1 00) pseudopure state 
using Eq.[3] and applied the quantum circuit of figure 7(a), for \x) = 1 00) , |01), 1 10) and |11). 
Finally, the spectra of proton and carbon were recorded individually in two different experiments 
by selective (7r/2) pulses after a gradient. The complete tomographed density matrices in each 
case is given in figures 7(d), (f), (h) and (j). In each case, the searched state \x) was found to be 
with highest probability. Ideally in a two-qubit system, probability should exist only in the searched 
state, and there should be no coherences. Experiment ally however, other states were also found with 
low probability, and some coherences were found in the off-diagonal elements of the density matrix. 
These errors are mainly due to relaxation and imperfection of pulses caused by r.f. inhomogeneity. 
Imperfection of r.f. pulses can cause imperfect refocusing of dynamic phase. However, it was found 
that setting the duration of selective pulses to múltiples of (2/J) yielded better results. We have 
used 13.2ms (6/J) duration Gaussian shaped pulses. The maximum errors in the diagonal elements 
are within 10% and that in the off-diagonal elemenst are within 15%. 

V. CONCLUSION 

A technique of using non-adiabatic geomètric phase for quantum computing by NMR is demon- 
strated. The technique uses selective excitation of subsystems, and is easily scalable to higher qubit 
systems provided the spectrum is well resolved. Since the non-adiabatic geomètric phase does not 
depend on the details of the path traversed, it is insusceptible to certain errors yielding inherently 
fault-tolerant quantum computation [45,46]. The controlled geomètric phase gates were also used 
to implement DJ-algorithm and Grover's search algorithm in a two-qubit system. Implementation 
of fault-tolerant controlled phase gates using adiabatic geomètric phase demands that the evolu- 
tion should be 'adiabatic', which requires long experimental time. To avoid decoherence, use of 
non-adiabatic geomètric phase might be utile. 
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Figure captions 

Figure 1: (a) The transport of a selected subsystem of two states \r) and \s) through slice circuit 
[18]. (b) Corresponding pulse sequence. A and B are transition selective pulses incorporated into a 
Hahn-echo [18], where A= (vr)P Hs> and B= (vr)^]. The path of the polarization vector under 
applied r.f. pulses is shown as A and B in (a). Due to the A= (ir)^^^ pulse, the polarization 
traverses a path from +z to -z. It comes back to +z along a different path if it is rotated about an 
axis with azimuthal angle (0 + n + 0), therby enclosing a sliced circuit of sòlid angle 20. (c) The 
transport of a selected subsystem of two states \r) and \s) through a triangular circuit [18]. (d) Pulse 
sequence for implementation of the triangular circuit given in (c). A=(iï/2)% C= (</>)j^ <_> '^ and 
B= (7r/2)g|^l' 5 j. The polarization vector is flipped to xy-plane by (tt/2)^^ s \ rotated about z-axis 
by (0)i r ^'^ and brought back to the z-axis by a (7r/2)l r >~l s > rotation about (0 + vr - 0). The sòlid 
angle enclosed by the circuit is 0. 

Figure 2: Observation of non-adiabatic geomètric phase when the subsystem is taken through the 
slice circuit of figure l(a). (a) Equilibrium 13 C spectrum of 13 CHC1 3 . The three small lines are from 
the natural abundant 13 C signal from the solvent of CDC1 3 . (b) The 13 C coherence of the prepared 
|00) + 1 10) state. (c) The 13 C coherence of the state 1 00) +e^\10). Since the subsystem goes through 
a closed circuit, the coherence gains a phase of purely geomètric origin, the magnitude of which is 
proportional to the sòlid angle (20) of the circuit. It may be noted that since all dynamical phases 
are refocused, the solvent signal from CDC1 3 is always absorptive irrespective of 0. The solute 
( 13 CHC1 3 ) signal on the other hand gains a phase of 0. The resulting signal changes shape from pure 
absorptive for = 0, to intermediate phase for arbitrary 0, dispersive for = 7r/2 and to absorptive 
(negative sign) for = ir. 

Figure 3: The geomètric phase gained by the 13 C coherence in Fig 2(c) is plotted against the 
corresponding sòlid angle (20). The observed phase closely matches the expected. 

Figure 4: Observation of non-adiabatic geomètric phase when the subsystem is taken through the 
triangular circuit of figure l(c). (a) Equilibrium 13 C spectrum of 13 CHC1 3 . (b) The 13 C coherence 
of the prepared |00) + 1 10) state. (c) The 13 C coherence of the state |00) + e"^ /2 |10), after the 
subsystem goes through the triangular closed circuit. The coherence changes from absorptive to 
dispersive and then to a negative absorptive with the change of 0. It may be noted that sign of 



phase of the coherence is opposite to that of 2(c). 



Figure 5: A plot of the absolute value of observed geomètric phase gained by the 13 C coherence in 
Fig 4(c) is plotted against the corresponding sòlid angle (0). The plot demonstrates high fidelity of 
such gates. 

Figure 6: Implementation of DJ-algorithm using non-adiabatic geomètric phase in the two-qubit 
system of 13 CHC1 3 . (a) Quantum circuit for implementation of DJ-algorithm in a two-qubit system 
[36]. h=(7r/2) w and h~ 1 = (n/2)_ y (b) Equilibrium l· l·l and 13 C spectra. (c) The l· l·l and 13 C spectra 
obtained after completion of the quantum circuit of (a) for Uf 00 , and application of a gradient pulse 
followed by (tt/2) pulses on 1 H and 13 C individually. (d) The complete tomographed density matrix 
[40]. The real and imaginary parts of the density matrix are given separately with the imaginary 
part being magnified five times (x5). (e), (g) and (i) are respective spectra obtained after C// u , 
Uf 0l and Uf 10 . (f), (h) and (j) are the corresponding tomographed density matrices. For constant 
cases (c) and (e), the final state is 1 00) , as shown in (d) and (f). For balanced cases (g) and (i), 
the final state in |10), as shown in (h) and (i). The diagonal elements have a fidelity of 95%, while 
off-diagonal parts have a fidelity of 87%. 

Figure 7: Implementation of Grover's search algorithm using non-adiabatic geomètric phase in the 
two-qubit system of 13 CHC1 3 . (a) Quantum circuit for implementation of Grover's search algorithm 
in a two-qubit system [42,30]. The U x and U 00 phase gates were implemented by non-adiabatic 
geomètric phase using selective excitation by 13.2 ms (Q/J) long Gaussian shaped pulses. (b) Equi- 
librium 1 H and 13 C spectra. (c) The l· l·l and 13 C spectra obtained after completion of the quantum 
circuit of (a) for \x) = |00), and application of a gradient pulse followed by (n/2) pulses on 1 H 
and 13 C individually. The intensities of the various lines in the spectrum gives the populations of 
the density matrix. (d) The complete tomographed density matrix [40] after implementation of the 
quantum circuit (a) for x = |00). The real and imaginary parts of the density matrix are given 
separately and the imaginary part is magnified five times (x5). (e), (g) and (i) are the spectra 
obtained when \x) = 1 01) , \x) = 1 10) and \x) = |11). In each case the searched state \x) was found 
with highest probability after implementation of the search algorithm with a fidelity more than 85%. 
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